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Figure Captions 



1. a) Mounting Arrangement on USNS Kane for NPS Equipment b) The Bouy De- 
ployed from USNS Kane 

2. Z* values for period 3/1 to 3/4 when KANE was anchored 80 miles west 
of Straits of Gibraltar. 

3. Z* values for period 3/10 to 3/12 when KANE was in survey area. 

4. Calculation of water vapor mixing ratio scaling parameter, q*, profile 
determination vs. bulk, bouy data only. 

5. Calculation of potential temperature scaling parameter, T*, profile de- 
termination vs. bulk, night time data. 

6. Calculation of friction velocity, U*, turbulence determination vs. bulk, 
all data. 

7. Calculation of evaporation duct height, Z*, profile determination vs. 
bulk. The open circles represent all data, the circled x's are the 
bouy data only. 

8. Duct height vs. relative humidity for typical meteorological condi- 
tions found during the Kane cruise*, a) sensitivity to wind speed 

b) sensitivity to sea surface temperature. 

9. Duct height vs. relative humidity for typical meteorological condi- 
tions found during the Kane cruise; sensitivity to the air-sea tem- 
perature difference (AT = T(Z) - T(0)). 
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ABSTRACT 

This is a report on atmospheric boundary layer measurements made 
by NPS personnel aboard the USNS Kane in February and March of 1978. 
The data has been analyzed to obtain Monin-Obukhov similarity para- 
meters and EM propagation properties using bulk, profile and turbu- 
lence data. The evaporation duct height, Z*, was found to vary from 
1.0 meters to 25 meters with a typical value of 16 meters. A compari 
son with profile and turbulence methods showed that the bulk method 
was accurate to within a factor of two, although some of the disagree 
ment has been attributed to profile measurement inaccuracies. A sen- 
sitivity analysis showed that the relative humidity and air-sea tem- 
perature difference are the most critical measurements for bulk calcu 
lations of Z*. 
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I . Summary 

A. Introduction 

This is a report on atmospheric measurements made by the Naval 
Postgraduate School Environmental Physics Group aboard the USNS Kane 
during operations in the North Atlantic in March of 1978. The primary 
goal of these measurements was to obtain measurements of the height of the 
radar evaporative (surface) duct from surface layer water vapor and 
temperature profiles and to evaluate the practicality of using a simpli- 
fied (the so called "bulk" method) method for more routine estimations. 

The bulk method requires measurements of only four quantities: the sea 

surface temperature and the air temperature, wind velocity and relative 
humidity at some reference height (usually 10 meters). It also has the 
advantage of considerably less stringent measurement accuracy require- 
ments. Due to the extreme difficulty in making accurate profile mea- 
surements at sea, the bulk method is essentially the standard for evap- 
orative duct height calculations. 

B. Description of the Problem 

In a dielectric medium, electromagnetic wave propagation is 
influenced by the refractive index. Just as an abrupt change in refrac- 
tive index causes EM waves to be "bent", so does a continuous gradient in 
refractive index. Since water vapor is a major contributor to N, there 
usually is a strong gradient in N over the sea surface because of evap- 
oration. This gradient causes EM waves to be refracted downward. If a 
horizontally propagating wave experiences a refraction greater than the 
earth's curvature, the wave is said to be ducted . Since the surface based 
duct is caused by evaporation, it is often called evaporation duct . The 
water vapor gradient (and therefore the M gradient) over the ocean typi- 
cally decreases with height in the surface layer. As a result, at some 
height, the N gradient has decreased below the value required for ducting, 
this is called the evaporation duct height , Z*. 

Since ducting can lead to signal enhancements of many orders of 
magnitude, knowledge of N and Z* in the surface is very useful. Fortunately, 
the physics of the surface layer is well understood and the specifications 
of the N profile and Z* can be reduced to the determination of only four 
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scaling paramters. This theoretical framework, known as Monin-Obukhov 
Similarity (MOS), is based upon the relation of the surface layer profiles 
of temperature, wind speed, water vapor and turbulence with the surface 
fluxes of momentum and sensible and latent heat. As one might expect, the 
problem is simple in principle but difficult in practice, primarily because 
the fluxes cannot be easily measured from an unstable platform. At present, 
the best technique is to measure the profiles at several different heights, 
solve the MOS expressions "backwards" to determine the scaling parameters, 
thus determining the N profile and Z*. This process, equivalent to infer- 
ring the surface fluxes from the profiles, is referred to as profile 
method . 

Multi-level measurements necessary for the profile method are 
not trivial. For example, from measurements at heights of 10 and 25 
meters above the sea surface, one can typically expect differences on the 
order of ^ 0 % in wind speed, 3 % in relative humidity, and 0.2C in tempera- 
ture. To determine the scaling parameters to a 10% accuracy, measurement 
accuracies required are: wind speed: 1%, relative humidity: 0.3% and 

temperature: 0.02C. Note that the distortions caused by the presence of 

the ship could exceed these accuracies. Atmospheric measurements even 
approaching this precision require highly trained scientific personnel 
with state-of-the-art instrumentation. These extremely taxing measurement 
accuracies have led to the interest in the bulk method. 

If one considers the MOS expressions for the gradient as mathe- 
matical derivatives, then the integral of the gradient from the surface 
(Z=Zq= 0) to some reference height (Z) specifies the variable (say tem- 
perature, T(Z)). Note that we have introduced two integration boundary 
conditions T(0) and Z^, where T(0) is the sea surface temperature and Z^ 
is the roughness length. If we assume that the wind speed at the surface 
is zero (U(0) = 0) and the relative humidity at the surface is 100% (H(0) 

= 100), only four quantities need to be measured: T(0), T(Z), U(Z), and 

H(Z). However, an unknown quantity Z^, which we have also introduced, 
must be empirically determined. 

Average temperatures measured during the cruise were 
T(0) = 16. 2C and T(Z) = 14. 7C, a difference of 1.5C - an order of mag- 
nitude greater than the two level example given previously. In short, the 



12 



bulk method has two advantages; fewer measurements are required and an 
order of magnitude greater toleration to error. This simplicity was 
obtained by introducing the quantity which has been empirically 
studied by statistically averaging vast numbers of over-water measurements 
of bulk, profile, and turbulence data. One of the primary questions about 
the bulk method is the applicability of an "average" to a specific 
situation. An even more serious question deals with the validity of 
extrapolating the MOS gradient expressions (essentially Z"^) down to the 
surface. 

C. Summary of Measurement Results 

The total data base for the cruise consists of 360 periods (half 
hour averages) of which 40 periods contained near surface (Z=2.8m) bouy 
profile measurements. A tabulation of the EM propagation data (Z*, Fm, 

N, dN/dZ and ) is given in Appendix C. There are several ways of 
looking at the comparison of the profile or turbulence values with the 
bulk values. For instance, the bouy profile values for the water vapor 
scaling parameter, q*, agreed very well with the bulk values, on the 
average, but the variance of individual values was about 50%. Since a 
considerable portion of this 50% is due to uncertainty in the profile 
measurements, these results are probably rather conservative. The com- 
parison is summarized in the following table: 

QUANTITY SCALING PARAMETER AVERAGE AGREEMENT RMS VARIANCE %VARIANCE 



Water Vapor 


q* 


excellent 


.12 gm/kg 


50% 


Temperature 


T* 


poor 


.045C 


50% 


Wind Speed 


u* 


fair 


.06 m/sec 


15% 


Duct Height 




fair 


10 M 


100% 



0. Conclusions 

Duct heights calculated on the basis of bulk measurements are 
probably within a factor of 2 of the correct value, and possibly much 
closer. The actual calculation of Z* from T(0), U(Z) and H(Z) can be 
done on a hand held programmable calculator. Given the sensitivity 
analysis of Section V, we feel that the real problem in routine appli- 
cations is the temperature measurements. The wind speed can be obtained 
with sufficient accuracy using "eyeball" estimates and a Beaufort scale. 
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The relative humidity measurements can probably be made with a sling 
psychometer from a well exposed portion of the ship. However, there is 
no escaping the fact that the air sea surface temperature difference 
must be measured to better than 0.5C. This requires well calibrated and 
exposed temperature sensors (with low drift characteristics). The air 
temperature measurement needs to be made with a good aspirator in a 
location reasonably free of ship influence. 

Typical bulk Z* values estimated during the cruise were around 
15 meters. Under high relative humidity conditions (rain and fog) Z* 
was as low as 1 meter. The maximum Z^ was around 25 meters. Consid- 
erable time variations were observed, mostly associated with frontal 
passages. In view of this, measurements should probably be made several 
times a day. 

One result of the theoretical analysis was the astounding 
volatility of Z* under stable conditions (air warmer than water). For 
the example with T(Z)-T(0)=1 .4C , a relative humidity of 90% gave Z* = 5 
meters whereas a relative humidity of 75% gave Z* = 60 meters. Unstable 
conditions showed much less variability. In fact, since the relative 
humidity over the ocean is rarely less than 40%, we expect Z* to almost 
always be less than 30 meters during unstable conditions. Therefore, 
stable conditions produce unusually strong surface ducting and the most 
variable ducting conditions. Although expected, the effect is larger 
than we would have anticipated. 

E. General Comments and Recommendations 

1. We feel another cruise is called for to clear up the 
uncertainties of the Kane data. We are now "on line" with improved 
temperature and humidity profile instruments, an IR sea surface tem- 
perature sensor and an improved computerized data logger and analyzer. 
One more cruise could do a great deal in establishing the bulk method 
procedures. 

2. An EM propagation analysis of other NFS data could be 
compared to the Kane data or to climatological predictions. We have on 
file data from the Pacific, the North Atlantic and the Mediterranean. 

3. More work needs to be done on the drag coefficients 
(related to Z^). The wind speed coefficient is fairly well known but 
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real knowledge about the temperature and water vapor coefficients is 
quite rare. 

4. Development of a simple system for routine shipboard 
measurements of sufficient quality to allow bulk calculations should 
probably be undertaken. 

5. A great deal of boundary layer work needs to be done 
(above the surface layer) to allow extension of similarity scaling to 
greater heights. If successful, it may be possible to make estimates of 
elevated ducts from surface and remote sensing data. 
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II. Instrumentation 



A. Sensors and Placement 

Rather than embark upon a detailed description of the NPS 
instrumentation, we will include as Appendix A a report (Houlihan et al , 
1977) on the subject and give only a brief sketch. The mean measure- 



ments are summarized in 


the following table: 




MEASUREMENT 


INSTRUMENT 


ACCURACY 


Temperature, T 


Quartz oscillator 


.OIC 


Relative Humidity, RH 


Li Cl cell 


• 3 % 


Wind Velocity, V=U 


Cup anemometer 


4% 



The relative humidity and temperature sensors are deployed in 

aspirated shelters to reduce errors from solar and IR radiation. These 

aspirators represent a weak point in the temperature profile measurements 

and do not produce reliable temperature profile data during the day. 

The turbulence data was based upon measurement of fluctuations of wind 

velocity (yielding the rate of dissipation of turbulent kinetic energy, 

e) and temperature fluctuations (yielding the temperature structure 

function parameter, Cj ). Constant temperature hot film sensors were 

used to measure e and microthermal platinum wires were used to measure 
2 

Cj . The frequency response of these measurements is greater than 
500HZ. 



The instruments are deployed at discrete levels: two levels 

on the NPS tower, one on the bow davit and (circumstances permitting) 
one on the buoy (see Figure la). The following are the level heights 
(Z) above the water: 



Level # 1 

Height, meters 2.8 



2 

9.1 



3 

18.6 



4 

24.7 



16 




Figure la. Mounting Arrangement on USNS Kane for NFS Equipment. 

The bouy, which is lowered into the water with a crane on the 
bow of the ship, is used to obtain a measurement very close to the 
surface. (See Figure lb). Since the gradients we wish to measure vary 
roughly as Z , they are easier to measure near the surface. Unfortu- 
nately, the buoy can only be used when a) the ship is stationary, 
b) the seas are moderate. 
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III. Theory 

A. EM Propagation 

1. Refractivity N and Ducting 

The refractivity N (in N units) of the atmosphere for EM 
wave propagation is 

N = 77.6 P/T + 3.73 x 10^ e/T^ (1) 

Where P is the atmospheric pressure, T the absolute temperature and e is 
the partial pressure of water vapor. The vertical gradient of N is 



M = r 
dZ '^1 8Z 



C + C — 
3Z ^3 3Z 



( 2 ) 



where q is the water vapor mixing ratio (q = .625 e) in gm/kg and 
C-| = .3, C 2 - 7.2 and - -1.3. A negative gradient in refractivity at 
height 7. causes a horizontally propagating EM wave to be bent towards the 
surface. If the gradient is large enough, the amount of bending will ex- 
ceen the curvature of the earth and the wave will be ducted at the height 
Z. Since the magnitude of the gradient decreases with increasing height 
in the surface layer, at some height (Z*) the gradient will no longer be 
strong enough for ducting. Waves propagating below Z* are ducted, those 
above are not. This height is called the evaporative duct height, and 
is defined by the critical gradient necessary for EM wave trapping. 



-,157 m-’ (3) 

Using Eq. 21 and noting that dP/dZ = -pg (the normal hydrostatic balance 
for the atmosphere), the Z* is defined by the following equation 

-.157 = -.032 + C^lf . C 3 fl (4) 

Since 3q/9Z and 3T/3Z are height dependent (approximately ''^Z"^) there 
will exist a value of Z = Z,* such that Eq. 4 is satisfied. 
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2. Structure Function Cj^ 

Turbulence in the atmosphere produces fluctuations in N on 
a broad range of size scales (milimeters to kilometers). These fluctua- 
tions scatter the EM waves producing a loss of resolution and coherence. 
If N(x) and N(x+d) represent values of N at two points, separated a dis- 
tance d, then 

= < (N(x) - N(x+d))^>d'^/^ (5) 

For EM Waves 



= ( 



t2 



) ' [C/ . I (1 . q ^ 



( 6 ) 



where A and B are constants is the water vapor structure function 

and C-j-q is the covariance function. Neglecting the 'Tq term (at sea level) 

. 2 _ r (77.6 X 10’®)(7.67 x 10^)-,2 
" L O J 



(7) 



T 



B. 



Monin-Obukhov Similarity Surface Layer Equations 
The vertical gradients of the wind speed (U), potential tem- 
perature (T) and water vapor mixing ratio (q) are given by 



^ (C) (8) 

Where X = U, T or q, k = Von Karmon's constant (35), is the diffusi- 
vity constant (a^j = 1.0), X* is the scaling parameter and (j)(c) is the 
MOS stability correction function (Appendix B). Also 

5 = Z/L (9) 

where L is the Monin-Obukhov stability length defined as 



L = 



^TU 

iTr; 



2 

.18 q*) 



(10) 



If we integrate Eq. 8 from Z = to Z, we obtain 

X(Z) = x(z„,) [in z/z^^ - U)] (11) 

Where 'l^y(c) is the profile stability function (Appendix B). 

2 

The structure function parameter C is related to X* by the 

A 

function f^(s) (Appendix B) 

Z’^/^ f^ U) (12) 

A similar form exists for e , the rate of dissipation of fluctuations 
of X 

. ><4- E,(c) (13) 

Where 6 is the Corrsin constant and E (c) is the MOS dissipation func- 

A A 

tion. In this format, = e. 
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C. MOS Parameter Determination Methods 
1. Profile Method 

The profile method is based upon Eq, 11, expressed in the 

following form 



X(Z) = ^ [In Z - t^(c)] + [X(0) 

X 



A 



(14) 



This can be expressed as a linear regression of the form 

X{Z) = mY + b (15) 

Where Y = In Z I'*' we have two or more levels of X(Z), we can fit 

the date to Eq. 16, the slope m, of the fit gives X* = a^icm. Note that 
it is not necessary to know the constant factor, b, (related to Z ) 

U A 

to find X*. The actual process is done iteratively in the following 
manner; 

a. Assume a value for L 

b. For each data point at height Z, calculate C and 

c. Least squares fit the multi-level data to obtain 

» 4"* and q^. 

9 

d. Calculate a new value of L from Eq. 10 

e. Return to b) until L converges 

Since the wind speed profiles from the Kane were consi- 
derably influenced by the ship, we used the bulk method to calculate 



U 



* ♦ 



2a Bulk Method 

The bulk method is also based on Eq. 15. If one knows 
values for Z , then X* can be calculated simply from X(Z) and X(0). 
Rewriting Eq. 15 



X* = 






Am - 

(In Z/Z 



X(0)) 



OX 



- t,JU) 



( 16 ) 



We see the form of the standard drag coefficient equation 

X* = (X(Z) - X(0)) 

Where is the drag coefficient for X 



^X 



4n 

[1 - (V)'' 



(17) 

(18) 
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In this form, is related to the neutral drag cofficient by 



'XN 



1/2 _ V 



<Z/Zox) 



( 19 ) 



The MOS stability parameter, can be obtained from 
Eq. 10 and Eq. 18. 



( 20 ) 



(1 - (ot-k) C 



-1 r 1/2 



T ’ ''TN 



K(?)) 



Where 



1/2 



5 =£5l 3N [(T(Z) - T(0)) + .18(q(Z) - q(0))] (2I) 

° ''' ^UN U(Z)^ 



In this case, the process is somewhat simpler, one cal- 
culates Cq from the data and solves Eq. 21 iteratively for c. Given c» 
the MOS parameters are calculated straightforwardly from Eq.s 18 and 19. 
3. Turbulence Method 

Given a reasonable estimate of 5 (from either of the two 

previous methods) one can obtain an independent value of the scaling 

2 2 

parameter from turbulence measurements of Cj , or e. For instance 
(Eq. 12) 

ij- = Z^/^/fj(?) (22) 

or (Eq. 13) , 

U/ = cZ/E|j(c) (23) 

D. EM Propagation Equations in MOS Form 

Near sea level, the gradient of M can be written (Eq.s 2 and 3) 



<{>t(5) 

dN/dZ = -.032 + (7.2 q* - 1.3 T*) 



(24) 



The duct height, Z*, is obtained from Eq. 25 by setting 
dN/dZ = -.157, the critical gradient for ducting. 

-(7.2 q* - 1.3 T*) 
ciyK (.125) 



Z* = 



(25) 



Given q^, T* and L, we solve Eq. 26 iteratively to obtain Z*. Note 
that no ducting will occur unless 

7.2 q* - 1.3 T* < 0 (26) 
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The refractive index structure parameter is calculated from Eq.s 7 
and 12 

,^2,^ (77.6X10-V.67x10^) j^ q.2 ( 27 ) 

IV, RESULTS 

A. Discussion 

Before presenting the results, it is worthwhile to consider 
several points about the absolute validity of the entire process. The 
various stability correction functions - 4>(C), f(C), E(5) - are empiri- 
cally determined functions. In fact, they are based on overland mea- 
surements. The NPS group has spent a great deal of effort in evaluating 
the validity of these functions over water and, so far, they appear 
to be quite good. There are other problems with some of the empirical 
constants and the drag coefficients. The profile method is dependent 
upon the product a^ic, both measured empirically. The estimates we 
have used are 0 |. = 1.35 and k = .35, from Businger et al (1971). The 
actual value of both constants is still a subject of controversy. 

The bulk method is dependent upon the drag coefficients 
*^UN’ ^TN ^qN‘ analysis, we assumed that water vapor has 

the same characteristics as temperature, that is 

“T ' "^q ’ ^OT " ^Oq ’ ^ *q^^^ 

There are theoretical reasons to believe that this is a good approxi- 
mation but not exactly correct. Under these assumptions, Cjj^ = Cqj^- 
For the actual form of the drag coefficients, we assumed Cyj^ had the 
wind speed dependence given by J. Kondo (1975). Based upon measure- 

ments summarized by Liu (1977) and NPS measurements of C.,. , we assigned 

-3 

a value of = 1.1 x 10 and solved for C as a function of Cq. 

Thus, for each data value of we directly calculated c from the for- 
mula 

5 = + •145^-'') C^<0 (23a) 
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B. EM Propagation Results 

A listing of the bulk calculations of the EM data is given in 
Appendix C (Table I-C). The data includes calculations of the Z=10 
meter values of N, dN/dZ and C|^ as well as Z* and Fm. Fm is the mini- 
mum frequency trapped by the duct as is calculated from 

Fm = 3.6 X 10^ (29) 

Where Fm is in GHz. 

The following are two examples to illustrate the general be- 
havior of Z*. Figure 2 shows Z* for a three day period at anchor ap- 
proximately 80 miles west of the Straits of Gibraltar. The steady de- 
crease of Z* for this first part of the record was associated with high 
relative humidity and rain. The rapid increase of Z* around 0200 hours 
on 3/3 corresponded to the passage of the front as indicated by clear- 
ing skies and a shift of the wind from 210^ to 260^. Figure 3 shows 
a very similar event on station in the survey area. In this case, the 
weather system was not strong enough to produce rain but fog was ob- 
served at 1800 hours on 3/10. The wind shifted from 200° to 300*^ be- 
tween 1900 hours on 3/10 and 0700 hours on 3/11. 




Figure 2. Z* values for period 3/1 to 3/4 when KANE was anchored 80 miles 
west of Straits of Gibraltar. 
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Figure 3. Z* values for period 3/10 to 3/12 when KANE was in survey area. 

C. Evaluation of the Bulk Method 
1. Scaling Parameters 

In order to make this evaluation more meaningful, we have 
analyzed the scaling parameters under those conditions of maximum profile 
accuracy. For instance, the humidity measurement were not accurate enough 
to give reliable profiles except when the bouy was deployed. The tempera- 
ture profiles were only analyzed for the night time periods. The turbu- 
lence measurements of U* were valid during all periods. A complete tabu- 
lation of comparisons for all periods is given in Table II-C (Appendix C) 
but keep in mind that the graphs shown below are based on data selected 
from the table. Incidently, the bouy data comparison has been extracted 
for condensed viewing and is shown in Table III-C (Appendix C). 

The comparison of bulk and profile calculations of q* are 
shown in Figure 4 from the bouy data. On average, the comparison is 
very good. The individual point scatter is about 50%. For this graph, 
and those following, the error bars are for the mean estimates and the 
number represents the number of data points in each bin. 
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Figure 4. Calculation of water vapor mixing ratio scaling parameter, 
q*, profile determination vs. bulk, bouy data only. 

The comparison of bulk and profile calculations of T* are shown 
in Figure 5. These results are very poor indeed. It is not clear if 
the lack of correlation is due to measurement problems or is in fact 
failure of the bulk method. 
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Figure 5. Calculation of potential temperature scaling parameter, T*, 
profile determination vs. bulk, night time data. 

The comparison of the bulk and turbulence calculations of U* are 
shown in Figure 6. Although the correlation is not perfect, it is pro- 
bably good enough for bulk determinations of Z*. The complicated velo- 
city dependence of that was used for these U* bulk values did not 
give significantly better comparison than assuming a velocity indepen- 
dent value of Cyj^ = 1 .3 x 10~^. 
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Figure 6 
tion vs. 



bulk and 
is given 




Calculation of friction velocity, U*, turbulence determina- 
bulk, all data. 

2. Duct Height 

Evaporation duct height, Z*, was calculated using the 
profile methods for all data. A comparison of these results 
in Figure 7. 
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Figure 7. Calculation of evaporation duct height, Z , profile deter- 

A 

ruination vs. bulk. The open circles represent all data, the circled 
x's are the bouy data only. 

Although the correlation is only fair, it should be noted 
that the bins containing most of the points (187 out of 232) fall 
within the expected mean error. Also encouraging is the much better 
agreement obtained from the bouy data. 

V. SENSITIVITY ANALYSIS 
A. Data Distribution 

Rather than consider many possible combinations of the me- 
teorological parameters, the sensitivity analysis will be done with- 
in the framework of typical values measured during the cruise. A 
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summary of the overall average data is given below. 



QUANTITY 


SYMBOL 


AVERAGE 


a_ 


Sea Surface Temp. 


TS 


16.25° C 


0.40° C 


Air Temperature 


T(10) 


14.67° C 


O 

o 

• 

1 — 


Relative Humidity 


H(10) 


70% 


17% 


Wind Speed 


U(10) 


7.57 m/sec 


2,84 m/sec 


Mixing Ratio 


q(10) 


7.29 gm/kg 


1.7 gm/kg 



B. Sensitivity of Z* 

The sensitivity analysis was performed using the meteorologi 
cal parameters as inputs into the bulk expressions. The predictions 
are shown in the form of Z* as a function of relative humidity. Of 
course, one should keep in mind the extreme rarity of relative humidi- 
ties (at Z=10 meters) less than 40% over the ocean. On each of the 
graphs shown, the solid line represents the "average" curve (TS-16.2, 
T(10)-14.7 and U(10)=7.6 m/sec). Since the average relative humidity 
was 70%, we expect an average Z* = 16 meters. Figure 8 shows the ef- 
fects of velocity variations on the order of the standard deviation 
found during the cruise. The effect of sea surface variation (at con- 
stant air-sea temperature difference) is given in Figure 8b. Z* is 
moderately insentive to both parameters. The corresponding graph for 
sensitivity to the air-sea temperature difference is noticeably differ 
ent (Figure 9). Not only can very large values of Z* occur for stable 
conditions (AT>0) but note the very steep variability of Z* between 
60% and 90% humidity. Of course, one must bear in mind that these are 
surface layer extrapolations, the maximum height of validity is proba- 
bly less than 50 meters. Nevertheless, these results are quite sig- 
nificant. 
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Figure 8. Duct height vs. relative humidity for typical meteorological 
conditions found during the Kane cruise; a) sensitivity to wind speed 
b) sensitivity to sea surface temperature. 
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Figure 9. Duct height vs. relative humidity for typical meteorological 
conditions found during the Kane cruise: sensitivity to the air-sea 

temperature difference ,(<^T = T(Z) - T(0)). 
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The large difference in stable and unstable atmospheric 
evaporation ducting properties is a direct result of the behavior of 
the MOS gradient equations. Let us assume that very near the surface 
we have an N gradient great enough to permit ducting. Since the gra- 
dient very near the surfaces decreases with height as Z \ we expect 
the gradient to decrease with increasing height until Z=Z* and the cri- 
tical gradient is reached. The atmospheric stability enters into the 
problem through its influence on the height dependence of the gradients. 

Under unstable conditions, the convective mixing quickly removes the 

-3/2 

gradients as Z increases, leading to a Z gradient height dependence. 
Consequently, the critical gradient is reached at a relatively low Z. 
Under stable conditions, the bouyancy forces are retarding the turbu- 
lent mixing as Z increases, leading to gradients that approach a con- 
stant independent of height. Therefore, the gradient decreases more 
slowly under stable conditions and the critical gradient can be reached 
at much greater heights. 

VI. CONCLUSIONS 

A. Variability of Z* 

Typical values of Z* obtained during the cruise ranged from 
around 2.0 meters to 25 meters with a mean value of 15 meters. Trans- 
missions at 10 GHz would have been ducted approximately 70% of the 
time, 20 GHz approximately 85% of the time. The two factors having 
the greatest influence on Z* are the relative humidity and the air-sea 
temperature difference. In one instance the duct height changed from 
3.0 meters to 20 meters in about four hours following a frontal passage. 

B. Validity of the Bulk Method 

The profile and turbulence data from this cruise do not give 
a clear validation of the bulk method. Although the water vapor and 
wind speed scaling parameters were reasonably well correlated in the 
mean, there was considerable scatter (on the order of 50% and 15% res- 
pectively) for individual measurements. Comparisons for T* and Z* were 
definitely a disappointment. We feel that much of the discrepancy can 
be attributed to inaccurate temperature and humidity profiles. Since 
both of these systems have been upgraded winch the Kane cruise, another 
cruise should be made. Based upon comparison of turbulence T* (from Cj ) 
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with bulk T* from other NPS data, we feel that the bulk T* are much better 
than indicated by the Kane data. 

C. Implications of the Sensitivity Analysis 

The sensitivity analysis showed that the wind speed and sea sur- 
face temperature contributions to Z*, though significant, were not highly 
critical. The relative humidity and air-sea temperature difference were 
definitely critical factors. Given the standard measurement capabilities 
of these two quantities, it is obvious that the air-sea temperature dif- 
ference will be the primary roadblock to implementing a program of rou- 
tine measurement of Z* either from ship data or satellite data. Another 
result of the study was the large values of Z* possible under stable con- 
ditions (stable conditions predominate off the east coast of the U.S. and 
Canada in the summer). 
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ABSTRACT 

Programs of investigation of optical propagation, marine fog, and 
boundary layer properties in the marine environment have led to the 
development of a shipboard system for measuring atmospheric stability 
and turbulence. This research requires very accurate multilevel 
measurements of mean and fluctuating temperature, wind velocity, and 
humidity. This paper is a description of the system including the 
following aspects: 1) instrumentation, 2) calibration, 3) analysis 

techniques, and 4) examination of the data and evaluation of the 
system. 
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I. Introduction 



There is currently considerable interest in curoulent transport pro- 
cesses in the marine boundary layer. Two problems which are of parti- 
cular importance concern some turbulence properties and their influence 
on optical propagation, and the correlation between the transport of heat 
and moisture and the occurrence of marine fog. In the optical problem, 
turbulence causes beam wander and spreading and scintillation (twinkling) , 
which cause energy density loss, resulting in information loss in the case 
of communication systems. The occurrence of marine fog is strongly depen- 
dent on the moisture/heat ratio in the atmosphere and the turbulent trans- 
port of these quantities is of critical importance. The goal of present 
investigations of turbulence properties is to enable predictions of the 
relevant atmospheric properties on the basis of relatively simple obser- 
vations that are currently being used by weather centers. A sophisticated 
series of experiments is necessary at this stage of development because 
of the need to develop an adequate parameterization of the problem and 
because of the complications inherent in model development and verification. 

The Naval Postgraduate School research vessel R/V Acania is the 
primary platform for overwater measurements and it is the system which is 
installed on this ship which is described here. A four-level system has 
been developed to measure both mean and fluctuating quantities. This 
system is augmented by a visiometer, a sea surface temperature probe, an 
acoustic sounder, and at times by balloons, kytoons or kites to make ele- 
vated measurements. 

Since a ship is a non-stable platform there is always some question 
as to the validity of data obtained when turbulence parameters are being 



investigated. It is not possible to measure turbulence directly at low 
frequencies (<lHz) aboard ship because the ship motion interferes in 
this region. At higher frequencies valid data can be obtained, and the 
experimental and analysis techniques described herein exploit this fact. 
An additional problem with shipboard measurements is the large physical 
size of a ship and the disturbance of the local air flow by the ship’s 
presence. Fortunately, the R/V Acania is a fairly small and narrow ship 
and this effect is not large. Moreover, proper placement of sensors to 
take data only when the wind direction is favorable can minimize inter- 
ference effects on the data. In this paper evidence is presented to show 
that neither ship motion nor physical presence invalidate the data 
displayed . 



II . Shipboard Configuration 

The shipboard sensor mounting arrangement is shown in Figure 1. The 
sensors for mean and fluctuating wind, temperature and humidity were 
mounted on two towers, located approximately 2 feet and 16 feet behind 
the tip of the bow. The heights of the sensors above the sea surface were 
Level 1 4.2 meters 

Level 2 6.6 meters 

Level 3 7.6 meters 

Level 4 13.9 meters 



The following sensors were installed at the levels indicated on the masts: 



Mean Wind 
Mean Temperature 
Mean Humidity 
Wind Fluctuation 
Temperature Fluctuation 
Humidity Fluctuation 



Cup Anemometer 
Quartz Thermometer 
Li Cl Cell 
Hot Wire 
Cold Wire 
Lyman - a 



All Levels 
All Levels 
All Levels 
All Levels 
All Levels 
Levels 1 and 3 
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The heights and locations of the sensors were chosen in such a way 
to minimize shipboard influence and to measure over a wide enough height 
difference so that small gradients in the parameters would be detected. 

A subsequent improvement that allows data to be obtained much closer to 
the sea surface will be described at the end of this section. 

The acoustic sounder was mounted on the fantail of the ship since 
this is the area that is most noise free. The vibrations of the deck 
when the ship is underway produce enough noise that useful signals can 
only be obtained up to about 400 meters. 

The visiometer was mounted on top of the bridge on a 3 feet high 
support. This position placed the instrument in an uncontaminated air* 
flow when the relative wind was within 90^ of the bow. 

The sea surface temperature sensor was located approximately 4 feet 
from the starboard side of the ship and floats within one inch of the 
surface. The starboard side was chosen because the ship water exhausts 
were located on the port side. 

In order to study near surface effects and to take advantage of 
larger gradients in the near surface region, a buoy system (Figure 2) 
was developed for deployment off the bow of the Acania (Corbin, 1977) . 

A heavy weight was suspended by cable through a hollow mast in the buoy, 
thus allowing the buoy to be rigidly controlled in the horizontal dir* 
ection but to freely slide up and down with waves. 

Ill . Instrumentation 

A. Mean System 

A block diagram of the mean system is shoxm in Figure 3. The 
system consisted of two primary components, viz., the sensors and their 
associated electronics, and the data logging and readout. 



The data logger was a micro programmable integrated data acquisition 
system (MIDAS) which was developed at NFS specifically for acquiring 
shipboard meteorological data. MIDAS utilized an Intel Model 8008 
Central Processor to control the sampling, averaging and recording of 
mean meteorological data. All software programs were written in PL/M 
to facilitate the writing of self -documenting programs (Plunkett, 1976). 

The operator was interfaced with the system via teletype for full 
duplex input/output communications and program control. The operator 
exercised control over the sample start-time and the number of samples 
to be averaged before outputting. The operator could also alter the 
preset sample list by adding or deleting various sensors as they came 
on line or become inoperable. Once initiated, the system was fully 
automated to sample the tailored list of sensors every 30 seconds and 
periodically printed output values averaged over a selected interval 
from one to sixty minutes. 

Output values were printed on the teletype in columnized format 
with the time of print as a leader. The teletype contained a paper 
tape punch that could be activated by the operator to produce a data 
copy concurrent with the printout. A magnetic cassette tape recorder 
was integrated into the system as a third data output device. This 
cassette is to be interfaced to a Hewlett-Packard Model 9831 Portable 
Computer so that profile and gradient flux estimates can be performed 
onboard automatically, using the BASIC programming capability of the 
HP 9831. Presently data cards are punched from the paper tape output 
and processed on the IBM 360 System at NFS. 

Mean wind speeds were measured using Thornthwaite Cup Anemometer 
units featuring lightweight plastic sensors that ensured low wind speed 



response. The pulse outputs (one pulse per revolution) of the four 
anemometer units were fed directly to MIDAS where they were averaged 
and printed on the teletype at the conclusion of each preselected time 
period. Revolution counts were also displayed on Hewlett-Packard Model 
5221A Counters for quick-look and calibration checks. The counts were 
converted to wind speed using the calibration parameters supplied by 
Thornthwaite, which were checked periodically in a wind tunnel. 

The temperature trnasducers utilized were Hewlett-Packard Model 2833 
Oceanographic (quartz crystal) Sensors, with a resolution of O.Ol^C. 

The sensors and accompanying booster oscillators were sequentially multi- 
plexed through a Hewlett-Packard Model 2801 Quartz Thermometer Readout 
Unit for monitoring. 

The humidity sensors utilized were Hygrodynamics, Inc., Model 1818W 
Wide Range Sensors, which featured an integral thermistor probe with a 
resolution of O.l^C. These units were sequentially multiplexed through 
a Hygrodynamics Digital II Readout Unit for monitoring both relative 
humidity and temperature data. 

The temperature probes in the humidity sensor arrangements provided 
a check capability on the more sensitive quartz thermometer sensors 
colocated with the humidity sensors. While the resolution demanded for 
gradient measurements (O.Ol^C) could only be satisfied by the quartz 
thermometer system, the availability of back-up, bulk temperatures from 
the humidity units was useful on several occasions where conflicting 
profile trends were discerned in quartz temperature readings. 

Temperature and humidity sensor outputs at each of the four 
measurement levels were multiplexed to their respective readout units 
using an NPS developed level selector which operated under the timing 



control of the MIDAS system. The level selector fulfilled two essential 
requirements. First, it allowed one readout unit to serve the entire 
profile array, thus decreasing the capital investment in equipment. 

Secondly, the level select unit contained the power elements for each 
level in the array. Thus, when individual sensors were selectively 
connected to their respective readout unit, only signal levels were 
switched at each junction. Maintaining constant powering at each level 
rather than switching in power concomitant with signal level switching 
kept line transients to a minimum and insured proper time constant 
response for the temperature and humidity equipment. 

The quartz temperature sensors and oscillator circuit, and the 
humidity sensors at each level were contained within C.C. BREIDERT CO., 
Air-X-Houster Aspirated Weather Shelters. The mounting arrangement is 
shown in Figure 4. The shelters served two purposes, viz., protecting 
the sensors from the marine environment, and reducing radiation effects. 
Radiation effects are especially severe in this application since some 
of the sensors are above the water and some above the deck of the ship, 
which is heated by the sun and radiates strongly. Because of radiation 
from the deck the shelters were modified by placing a shield at the base 
of the unit to reduce radiation from below. 

B. Fluctuation System 

A block diagram of the fluctuation system is shown in Figure 5. 

It is immediately apparent that there are several output devices, allow- 
ing several methods of data analysis. The analog tape recorder was used to 
record the basic signals, which could then be analyzed in the laboratory 
subsequent to a field trip. The analog -digital converters and the 
averaging and storage were contained in the MIDAS system described above 



and their results were printed on the teletype at the same time as the 
mean system data. Averaging could also be accomplished directly on the 
strip chart records. 

All sensors were mounted on wind vanes, with one vane per level, 
to maintain correct alignment with the relative wind directions. 

The sensors for wind speed fluctuations were Thermo -Systems model 
1210W-T1.5 probes (4.5 micron tungsten) operated as hot wires with a 50% 
overheat. The power and bridge for the probes was the Thermo-Systems 
model 1054 constant temperature anemometer, the output of which was 
processed by a Thermo-Systems model 1057 signal conditioner to remove 
the dc level. The dc level could cause the tape recorder to saturate 
when large fluctuation signals occured. The signal conditioner also 
contained band pass filters, which were set for a range of 0.2 Hz to 
1 kHz, reducing electronic noise. The 50% overheat used was high enough 
so that the system response to temperature fluctuations was well below 
the electronic noise level. The hot wire sensors were oriented with the 
axis of the wire in the vertical direction. With this orientation the 
sensor was insensitive to the vertical component of the wind, so that 
only horizontal fluctuations were recorded. 

The sensors for air temperature fluctuations were Thermo-Systems 
model 1210W-P.8 (2.5 micron platinum) operated as cold wires. An ac 
resistance bridge, Sylvania model 140 thermosonde, which operated at 
3 kHz, was adapted for shipboard use and acted as the detector. With 
this unit data could be obtained up to 1 kHz. It is necessary to use 
an ac bridge as the detector so that high sensitivity could be obtained 
with a very low current through the sensor. If the temperature of the 
sensor was elevated due to bridge current the system would sense wind 



speed fluctuations as well as temperature f luctuat ions • The bridge had 
a long time constant, automatic balance circuit so that slow drifts of 
ambient temperature were compensated. 

The thermosonde configuration was such that external resistances 
would be used in both arms of the bridge. This allowed one to use a 
single sensor balanced by a precision decade resistance, or to use two 
cold wire sensors in the bridge. Thus, two methods were available for 
evaluating temperature fluctuations: (1) A single sensor, the signal 

being analyzed by a spectrum analyser-T* mode, and (2) two sensors spaced 
at a fixed separat ion-AT mode , the resulting signal being processed by 
an RMS module, with the average either recorded on a strip chart or by 
the MIDAS system. 

Figure 5 shows that the wind speed fluctuation signals were also 
processed by an RMS module after the signal had passed through a band 
pass filter. The band pass filter operated in the time domain and per- 
formed a function equivalent to the spatial filtering accomplished by 
the probe separation utilized in the AT mode for temperature fluctuations. 
Spatial filtering could not be used for wind speed fluctuations because 
it was not possible to operate two hot wires by a single bridge. Using 
two bridges and a different circuit introduced too much error into the 
process . 

A final arrangement for evaluating fluctuation signals led to the 
differentiation of the bridge output prior to FMS analysis. Dif- 
ferentiation necessarily emphasizes high frequency components so that 
high frequency noise becomes a problem. To alleviate this problem a 
1 kHz active filter with a 48 dB rolloff was used but noise problems 
were still severe enough so that no satisfactory results were obtained 
with this method. 
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Methods for analyzing data gathered utilizing the various experi- 
mental techniques will be described in Section VI. 

Humidity fluctuations were detected with a Lyman-a cell. Difficul- 
ties were experienced in keeping this system in operation in the rigors 
of a marine environment and no results will be reported here. 

C. Auxiliary Equipment 

An Aerovironment model 300 Acoustic Radar was installed on 
the R/V Acania in order to monitor the height of the temperature inver- 
sion. The sounder transmitted a 1600 Hz pulse of sound which was 
scattered by small scale ambient temperature inhomogeneities. Using a 
facsimile recording device, the scattered return sound waves were then 
detected to produce the observed signal. The large temperature gradient 
and strong turbulence at the base of an inversion produced a strong 
reflected signal, so that the height of the inversion could be readily 
determined. It was possible to measure the strength of thermal turbu- 
lence with the acoustic sounder, but noise and calibration problems made 
such measurements impossible with our installation. 

A Meteorology Research, Inc. model 1580A Fog Visiometer was installed 
on the bridge of the R/V Acania. The instrument measured the amount of 
light scattered from particulates and had a usable visibility range of 
8 m to 13 km. The unit could be rotated so that the relative wind was 
incident directly on the open side of the scattering region. Rotation 
was necessary when the sun was low on the horizon in order to keep sun- 
light out of the receiving optics. 
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IV. Calibration 



A. Mean System 

1. Temperature 

The mean temperature sensors (HP model 2833 Quartz) were cal- 
ibrated using standard laboratory thermometric techniques. Since the 
system was used to measure the gradient of a temperature profile, the 
relative consistency of the sensors was considered more important than 
the absolute temperature calibration. All five sensors were greased and 
screwed into a 3 inch aluminum cylinder fitted with a center -mounted 
thermocouple. The aluminum cylinder was 0-ring fitted into a copper 
pipe long enough (24 inches) to permit coiling of the majority of the 
armored sensor cable in the temperature controlled environment. The 
pipe and cylinder unit were immersed in a 15 liter glass dewar filled 
with a mixture of water and ice. The cylinder and pipe combination 
allowed the sensors and cables to be surrounded by the bath without 
actually getting wet. The large mass and excellent thermal conductivity 
of the copper and aluminum eliminated transient fluctuations and reduced 
inter-sensor temperature differences to a few thousandths of a centi- 
grade degree. 

The normal calibration procedure was to cool the sensors to about 
15°C while monitoring their indicated temperatures with printout from 
MIDAS at one minute intervals. Once equilibrium was reached, the sensor 
oscillators were adjusted so that (if possible) all sensors read within 
.01°C of each other. The temperature was then varied in 2-3°C incre- 
ments by adding ice or warm water and noting the MIDAS printout values 
when a new equilibrium was reached. Initially, the temperatures were 
compared to the thermocouple values (data set 1/7/76 in Table I) but 
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1/7/76 



2/9/77 



TABLE I Mean Temperature Calibrations 
T, ''C 



TC 


0 


1 


2 


3 


4 


.01 


.02 


.00 


.02 


-.01 


.02 


5.25 


4.94 


4.93 


4.94 


4.92 


4.94 


7.50 


7.24 


7.24 


7.25 


7.23 


7.24 


9.45 


9.21 


9.22 


9.22 


9.21 


9.21 


11.95 


11.64 


11.66 


11.65 


11.65 


11.65 


13.35 


13.05 


13.08 


13.06 


13.06 


13.06 


14.85 


14.53 


14.56 


14.54 ■ 


14.55 


14 .54 


23.25 


22.86 


22.89 


22.87 


22.88 


22.87 




8.61 


8.60 


8.62 


8.61 


8.62 




13.14 


13.14 


13.15 


13.15 


13.16 




20.54 


20.50 


20.51 


20.52 


20.53 


this 


was later deemed unnecessary. 


One problem with the system developed 


due to using quartz 


sensors with different 


calibration 


properties. In 


data 


set 1/7/76 of 


Table I, sensors 


1 and 3 


are from a 


different batch 



than 0, 2 and 4. This led to increasing disagreement at the extremes 
of the calibration range. This problem was alleviated by ordering more 
sensors and specifying characteristics to match sensors 0, 2 and 4 (data 
set 2/9/77 of Table I), With this improvement, the accuracy of a single 
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sensor (relative to the average of all five) was no worse than ± .02°C 
over a reasonable operating range. Note that this is the accuracy 
of the temperature sensor only and does not include errors introduced by 
the aspirators. The relative aspirator temperature error depends most 
strongly on the amount and uniformity of solar heating; this error is 
quite small at night or under heavy clouds but may be as great as 
± .l^C in bright sunlight with light winds. 

2. Humidity 

To determine their accuracy and reproducibility, detailed exam- 
ination in laboratory conditions has been made of the calibration char- 
acteristics of the Hygrodynamics Digital Hygrosensor. The calibrations 
were based on ambient humidity variations (from 24% RH to 50% RH over a 
one month period) , as measured with a Bendix Psychron Wet-Dry Bulb 
Psychrometer , and controlled humidity (from 40% RH to 99% RH) , as meas- 
ured by a wet -dry thermocouple ina2x3x3ft fog chamber adapted for 
these experiments. The air circulated through the chamber could be pre- 
humidified in a column of wet glass beads giving nominal humidity of 
about 80%. High humidities (up to 99% RH) could be obtained by intro- 
ducing varying amounts of water in the form of a fine mist into the 
chamber. In general, the humidity in the chamber could be stabilized 
for periods of several hours and was uniform throughout the chamber to 
within 1% RH. 

The hygrosensor was a Dunmore-type lithium chloride sensor whose 
resistance varies in proportion to the relative humidity to which it 
is exposed. It is a slow-response sensor and is generally utilized in 
mean system measurements. The standard procedure of calibrating these 
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sensors in a small chamber with a saturated salt solution of known vapor 
pressure was not .satisfactory for a program of calibration of eight sen- 
sors. This was primarily due to a lack of confidence in our ability to 
reproduce humidity conditions in consecutive sensor calibrations. Attempts 
to calibrate all eight sensors simultaneously in a large volume chamber 
were complicated by temperature drifts that made it difficult to assure 
equilibrium conditions necessary to the accuracy of the technique. As 
a result, it was decided to use the humidity controlled chamber previously 
described. The quantitative results of a series of calibrations are given 
in Table II; in general, the eight-sensor average agreed with the psychro- 
meter standard to -0.4 ± 2.9% RH for measurements taken with equilibrium 
times of order one hour. The individual total humidity accuracy of the 
sensor was of less importance than the relative consistency within the 
group of sensors, since they were used in a four-level system to determine 
humidity gradients. In this respect, they could be calibrated in the 
laboratory to ± 1% RH, subject to assurances of sufficient time for 
response. 

The suitability of the hygrosensors for shipboard measurement of 
humidity gradients was reduced by two inherent properties of the device. 

An aspirated, teflon-coated sensor has a time constant of about 30 minutes, 
but this varies from sensor to sensor. Therefore, the humidity differ- 
ence read between two sensors will be unreliable if the humidity is 
changing faster than about 10% RH per hour. A more serious deficiency 
appears if the sensors are exposed to humidities above 95% RH. At those 
high humidities, the sensors become sluggish and exposure to 99% RH may 
cause the sensors to be useless for as long as several days. Further 
uncertainties are introduced by the use of long cables in the shipboard 
system and by the accumulation of sea salt on or near the sensors. 
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TABLE II Hygrodynamics Digital II Calibrations 



Wet -Dry Ave 

Bulb Hygro 

%RH %RH 


1 


2 


3 


4 


5 


6 


7 


8 


33. 


31.6 


1.4 


-1.4 


-0.6 


1.5 


-0.2 


-1.6 


-0.9 


1.8 


36. 


37.2 


-0.5 


-0.8 


-0.3 


1.6 


0.5 


-1.4 


-0.1 


1.2 


44. • 


41.2 


-0.8 


-0.4 


-0.3 


0.6 


0.9 


0.4 


-0.1 


0.3 


NA 


41.5 


-0.5 


-0.4 


-0.3 


0.5 


0.9 


-0.4 


-0.1 


-0.1 


NA 


42.6 


-1.2 


-0.5 


-0.5 


-0.9 


1.8 


-0.4 


-0.2 


0.1 


NA 


65.8 


-0.2 


-1.4 


-2.5 


0.8 


1.8 


-0.3 


1.3 


0.3 


74.8 


73.6 


0.2 


0.2 


-1.3 


0.5 


-0.2 


0.4 


-0.4 


0.8 


88.3 


85.0 


-0.7 


0.9 


-0.4 


0.1 


-0.6 


0.0 


0.4 


0.2 


82.6 


86.7 


-1.4 


-0.6 


-0.5 


0.1 


-0.2 


0.3 


0.7 


0.3 


82.6 


87.0 


-1.0 


0.5 


-0.7 


0.0 


0.3 


0.7 


0.0 


0.2 


85.0 


83.0 


-0.4 


1.1 


-0.1 


0.1 


-1.8 


1.1 


-0.6 


-0.1 


93.9 


89.6 


0.1 


0.8 


-0.2 


0.3 


-2.2 


2.1 


-1.5 


0.4 


94.4 


92.3 


-0.6 


1.1 


-0.5 


-0.3 


-0.4 


1.2 


-0.5 


0.4 


94.4 


93.6 


-0.8 


1.0 


-0.8 


-0.3 


0.4 


-0.2 


0.8 


-0.4 


Ave 


Deviat. 


-0.5 


-0.0 


-0.6 


0.3 


0.1 


0.1 


-0.1 


0.4 




O 


0.7 


0.9 


0.6 


0.7 


1.1 


1.0 


0.7 


0,6 


The average 


deviation 


of the 


wet -dry 


bulb from 


the 


Ave Hygro 


reading 





0,4% RH with a standard deviation of 2,9% RH. 
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3. Velocity 

The mean wind velocity sensors (Thornthwaite cup anemometers) 
were calibrated in a wind tunnel with a nominal one meter square test 
section. The cups were calibrated four at a time with a symmetrical 
placement to reduce relative velocity differences. A summary of the 
calibration is shown in Figure 6 expressed as individual cup velocity 
deviations from the wind tunnel velocity. From this data it can be 
concluded that, for gradient purposes, the cup anemometers were accurate 
to about ± 4%. The deviation of the calibration curve at low wind 
speeds is due to errors in the determination of wind tunnel velocity 
values . 

B. Fluctuation System 

The essence of the fluctuation calibrations is contained in 
establishing the relation between the instrumental voltage variation, V*, 
about mean voltage, V, with respect to the atmospheric fluctuation, X*, 
about its mean value X, i.e. 



X* = F ( V, X ) V’ 



( 1 ) 



The calibration factor F ( V, X ) may be a function of individual 
sensors as well as the instrument being used. This process is usually 
accomplished by measuring the dependence of X on V and relating it to 
F ( V, X ) through the derivative. 



F ( V, X ) 




( 2 ) 
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1. Temperature 

The appropriate temperature fluctuation expression is: 

T' = ( a G R )”^ V (3) 

where a is the temperature-resistivity coefficient of the platinum sensor 
wire, G is the gain factor of the thermosonde (set at G = 20 volt/H) , and 
R is the resistance of the sensor. In a previous work (Schacher and 
Fairall, 1976) it was found that a = .0036 ^ with little variation 

from sensor to sensor. The frequency response of the thermosonde units 
was about 1 kHz. 

2. Humidity 

The Lyman-alpha sensor measures humidity as a function of the 

0 

absorption of ultraviolet light at the 1215A Lyman-alpha transition of 
hydrogen in water vapor. It consists simply of a UV source tube and 
detector, separated by an absorbing air gap. 

The detector voltage is proportional to the total light trans- 
mitted and is related to the absolute humidity, q (in mb), by the 
equation 

V = exp (-Xq) (4) 



The Lyman-alpha is a fast response device and is specifically used for 
measurement of humidity fluctuations, q’ 



q = 




11 

V 



(5) 



The quantity in Eq. (4) is the dry-air voltage which can be found by 
bathing the detector in dehumidified air. Since water soluble windows 
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(LiF or MgF^) are required for transmission in the ultraviolet, is a 



highly variable quantity, a property which makes the Lyman-alpha a poor 
device for measurement of total humidity, q. However, Eq. (5) shows that 
the fluctuating part, q’, depends only on X, which is a property of water 
vapor and the gap setting only and is independent of the window condition 
so that q’ can be reliably measured. The Lyman-alpha was calibrated in 
dry air, ambient air, and in the humidity chamber and it was found that 
X = .225 mb ^ for source = 50 ]ia and a 1 cm air gap (Fig. 7). 

Feasibility of the Lyman-alpha sensor for shipboard humidity fluc- 
tuation measurements is subject to considerations of the survivability 
of the windows in the ocean environment of high humidity, sea spray, and 
rain. In the laboratory humidity chamber, was found to be reduced by 
a factor of four after a three-day exposure to 85% humidity. Consequently, 
it was necessary to provide physical protection for the sensors during 
routine operations. 



calibrator. It was important to use the actual cables from the shipboard 
system (or their equivalent) to prevent errors due to voltage drops in 
the cable. The appropriate mean dependence function is 



3. Velocity 



The hot wires (or hot films) were calibrated using a TSI 1125 



V (U) = (V 



o 



2 




( 6 ) 



2 

where V and B are obtained for each sensor at a specific overheat. The 



o 



derivative of (6) yields the fluctuation calibration factor 



U 




(7) 



B 
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The frequency response of the velocity bridge was about 20 kHz. The 
sensor frequency response depends on the sensor diameter and overheat, 
but was typically greater than 1 kHz (Fairall and Schacher, 1977). 
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V. Data Acquisition 



Three criteria applied when deciding the methods for acquiring data: 

1) data was to be taken for a wide range of local meteorological con- 
ditions, 2) a given condition should be tested several times, preferably 
on different days and/or locations to test the statistical validity of the 
data, and, most importantly, 3) shipboard conditions must be such that 
valid data can be obtained. 

The most important consideration to insure validity of the data was 
the airflow over the ship. Since all of the sensors were located close 
to the bow of the ship the relative wind had to be from the bow or the 
ship could severely influence the data obtained. It was found that 
data could be taken when the wind was within 30° of the bow. When the wind 
was more than 45° from alignment with the bow the wind profiles were 
noticeably affected. The position of the sun was a second consideration. 
When the sun was low on the horizon and behind the ship the aspirators 
for levels 1 and 2 were in shadow and those for levels 3 and 4 were 
in sunlight. Under this condition the temperatures at the upper levels 
were elevated enough to make the obtained temperature profiles invalid. 

Since data was acquired on a fluctuating system, a time average had 
to be performed to obtain final parameters. Averaging over 10 minute 
and 20 minute time spans was accomplished depending on the situation. 

Thus, data had to be acquired in such a way that the conditions remained 
constant over a time at least as long as the planned average time. This 
means that there had to be close coordination between scientific personnel 
and ship crew so that no course changes occurred during a data taking run. 
Of course, changes in local conditions also affected the data, such as 
moving through the edge of a fog bank. Even seemingly inconsequential 
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occurrences such as moving into the lee of an island when the ship was 
several miles offshore could affect the data. A running log of all 
ship maneuvers and position, and of local conditions was maintained and 
no averages were performed for which it was suspected that any condition 
changed during the time period under consideration. 

Several tests were made and it was found that data could be taken 
both with the ship at rest and underway. However, some of the subsequent 
calculations then required that the true wind over the surface of the 
water be known accurately. Uncertainties as to the ship’s actual speed 
make calculation of true wind from the relative wind and ship heading 
and speed too inaccurate. Thus the ship was stopped frequently in order 
to obtain the true wind. 
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VI. Data Analysis 



The analysis procedures changed considerably as the instrumentation 
evolved towards increased automation. Averaging of mean data, initially, 
a painstaking hand process, is now performed by MIDAS at specified time 
intervals (typically 10, 20 or 30 minutes) . Rather than describing all 
techniques that have been used, discussion is limited to those of more 
recent application. 

A. Mean Profile Analysis 

The mean profiles and gradients were established under the assumption 
of the logarithmic height dependence. 



X = A log Z + B 



( 8 ) 



The gradient at the height Z is 



= A 

dz Z 



(9) 



In the case of temperature, one is interested in the potential temperature, 

0 , 

■0 = T + .0098 Z (10) 

and the virtual potential temperature, 6^, 

W = J + .61 q f (11) 

V 

where q is the specific humidity in grams of water vapor per gram of dry 
air. The actual fits of the data were done either by the standard least 
squares method or by subjective ("eyeball”) graphical techniques. The 
presence of individual erroneous values due to printer errors, micropro-^ 
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cessor errors, and instrumental glitches required the data to be subjec- 
tively edited even in the least squares analysis. In either case, one is 
required to identify and reject ”bad” values - a delicate and unsatisfying 
process. Figure 8 shows examples of profiles in temperature, humidity and 
velocity. 

B. Turbulence Analysis 

There are four different methods available to calculate turbulence 
parameters from the fluctuations in the atmospheric variables: spectral, 

derivative (dissipation), difference, and RMS. The spectral method is 
based upon the assumption of ”local isotropy” and the Kolmokorov -5/3 
slope of the one-dimensional power spectral density, , 

7 - 5/3 

4) (k) = .25 C (12) 

X X 



where is the structure function for the parameter x and k is the wave 
number. Using Taylor *s ”frozen turbulence” hypothesis (k = Zirf/U), C 

X 

is found by performing a fourier spectrum analysis of a signal in the fre- 
quency domain (f ) , 



f (f) 



k 4) (k) = .25 C ^ 

X X 




(13) 



Therefore, 

(14) 

which is expected to apply in the frequency range from about 0.1 to 100 Hz. 
(This is the inertial subrange of turbulence.) The structure function can 
also be found by measuring the variance of the difference in the variable 
X at two points separated by a known distance, d, 

C ^ = [x (r) - X (r + d)]^ (15) 

X 
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is independent of d in the inertial subrange. 



Defined in this manner, C 



X 



The RMS method is based upon measuring the variance of the signal fluc- 



tuations between selected frequency limits and relating this to through 



Eq. (13): 



k 
/ ^ 



4> (k) d k = X'^ = (X' 

X rms 



(16) 






Again assuming k = 27rf/U, then 



8/27T12/3 






(17) 



where f^ and f^ are the upper and lower frequency limits, and X* 

represents the measured RMS fluctuations of the filtered signal. Since 

the rate of velocity turbulence dissipation, £, is traditionally used, 

2 

rather than C , one should use the relation 
u 

= 2.0 (13) 



to calculate £• Due to the extreme sensitivity to high frequency noise, it 

was found that the derivative method was unsatisfactory for shipboard data. 

2 2 

The spectral method was used to calculate , C and £ (data shown 

1 q 

in Figure 9) . This method has the advantage that the quality of the data 
is immediately obvious and noise spikes (60, 120, 130 Hz) can be ignored. 

The primary disadvantage of the spectral method is the time required to 
analyze four levels of data with two or three signals for each level. Al- 
though it requires two microthermal sensors per level, the difference method 

2 

has proven to be an ideal technique for measuring because the RMS temp- 
erature difference fluctuations (Eq.l5) can be automatically averaged using 
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the ADC capability of MIDAS. Since this is a temperature difference 
measurement, it is less subject to external sources or error (such as 
humitity - salt induced temperature spikes (Friehe, 1977)). The RMS 
technique has proven to be quite satisfactory for measuring £, since it 

too lends itself to automatic averaging on MIDAS. Comparison of spectral 

• 1/3 

and RMS measurements of £, through the quantity = (K Z £) , 

(Figure 10) are quite consistent. is the friction velocity, K is 

Von Karmon’s constant (.35) and Z is the height above the surface. This 

expression is valid for near neutral stability where most of the data 

was taken. 
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VII . Experiments Performed 



Observational experiments were conducted aboard the R/V Acania 
(Fig. 1) during cruises to San Nicholas Island (SNI) , in the vicinity 
of SNI, in Monterey Bay, and on an open ocean cruise in the fall of 1976. 
Dates and locations of the experiments which yielded results presented 
here are listed in Table III. 

The route taken to and from SNI appears in Figure 11 for that portion 
of the cruise for which data was taken. The ship’s anchorage, position, 
and track in the vicinity of the northwest tip of SNI appear in Figure 12. 
For Monterey Bay results, all data were collected while the R/V Acania was 
anchored. The general anchorage locations which were used in Monterey 
Bay appear in Figure 13. The anchorage locations were changed to accom- 
modate coincident ship to shore or shore to shore optical propagation 
experiments. Wind directions during the Monterey Bay experiments were 
generally from the north to northwest so the location of anchorage is not 
expected to be a factor in the results and they are not distinguished with 
regard to location. Details pertaining to the 1976 open ocean cruise appear 
in a separate report, Fairall, et al. (1978). 

Changes were made in the shipboard measurement systems and sensor 
mounting arrangements during the sequence of experiments associated with 
these results. In general, mounting arrangements were changed to improve 
the vertical resolution of the measured parameters. During the first SNI 
cruise (1-4 March 1973), fluctuation (temperature and velocity) measure- 
ments were made at only two levels and mean wind, temperature and humidity 
measurements were made only at one level. Sea surface temperature was also 
measured. During the second SNI cruise (19-23 September 1973), both fluc- 
tuation and mean measurements were made at four levels which differed 



A-27 



slightly from the final configuration shown in Figure 1. All the Monterey 
Bay experiments were performed with the final mounting configurations 
described above. Sensors described for fluctuation and mean measurements 
were used in all experiments but the recording systems, particularly that 
for the mean data, were changed. 

Additional discussions on analyses, instrumentation and experimental 
conditions, along with additional results, are contained in theses listed 
in the references, which utilized data from these experiments. 
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TABLE III 



SUMMARY OF OBSERVATIONAL EXPERIMENTS TO FEBRUARY 1977 





Date 


Location 


Prevailing Condition 


1973 


1-5 Mar 


SNI 


Unstable 




18-21 Sep 


SNI 


Neutral-Stable 




15-19 Oct 


Mtry Bay 


Unstable 




26 Nov-5 Dec 


Mtry Bay 


Unstable 


1974 


15-18 Jan 


Mtry Bay 


Unstable 




25-26 Feb 


Open Ocean 


Stable 




27-29 Feb 


Mtry Bay 


Stable 




6-7 Mar 


Mtry Bay 


Stable 




13-15 Mar 


Mtry Bay 


Stable 




25-28 Mar 


Mtry Bay 


Unstable 




17-21 Jun 


SNI 


Neutral 




12-16 Aug 


Mtry Bay 


Neutral-Unstable 




16-21 Sep 


Open Ocean 


Unstable 




18-22 Nov 


Mtry Bay 


Neutral-Uns table 


1975 


9-10 Jan 


Mtry Bay 


Stable 




24-28 Mar 


Mtry Bay 


Unstable 




20-24 May 


Mtry Bay 


Neutral 




25-29 Aug 


Mtry Bay 


Neutral 




11-16 Sep 


Open Ocean 


Unstable 




9-11 Dec 


Mtry Bay 


Unstable 


1976 


29 Mar-2 Apr 


Mtry Bay 


Unstable 




26-30 Apr 


Mtry Bay 


Neutral-Unstable 




26-30 Jul 


Mtry Bay 


Neutral 




20 Sep-14 Oct 


Open Ocean 


Neutral-Unstable 


1977 


15-18 Feb 


Mtry Bay 


Neutral 
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VIII, System Evaluation 



A. Mean System 

Accurate measurements of atmospheric surface gradients, especially 
in the marine boundary layer, are extremely difficult. The severity of 
the problem can be illustrated by noting typical conditions in the marine 
boundary layer (Table IV) in terms of the typical gradients (8X/9Z at 
Z = 10 meters), scaling parameters (X^) and total difference from Z = 5 to 
Z = 15 meters (AX = |x (15) - X (5) [). The scaling parameter is defined by 

where f is a stability corrections factor (Businger, et al., 1971). To 
measure X^ to 10% accuracy requires the single level accuracies indicated 
in the table, compared to our present accuracies. The temperature, T, is 
used in the table rather than the virtual temperature because it is the 
variable actually measured. 

Although the gradients are not required to calculate the turbulent 

parameters, they are of utmost importance in evaluating the theoretical 

behavior of the boundary layer through specification of the stability 

(Richardson number) and through normalization of the turbulence parameter, 

2 

(Wyngaard, et al. , 1971). This lack of gradient accuracy led to con- 
struction of a near surface buoy system to exploit the Z ^ dependence of 
the gradients. With the present systems, the gradient measurements repre- 
sent the weak link in efforts to evaluate the turbulent boundary layer 
similarity expressions in the marine environment. 

B. Fluctuation System 

Despite a large array of factors influencing the performance of the 
turbulence measurement system (noise, ship effects, salt on sensors. 
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TABLE IV 



Typical Marine Boundary Layer Profile Parameters 



T, 



U/ m/sec 



q f gm/kgm 



Ox/32) 



10 met 



- .020 



.040 



- .025 



- .07 



.15 



.08 



Ax 



.20 



10 % 



5 % 



Present Accuracy 
Required Accuracy 



.04 



5 % 



2 % 



3 % 



1 % 



A- 31 



calibration accuracy) the values of turbulence parameters were quite con- 

2 

sistent. Figure 14 is a time history of four levels of and z showing 
the expected height dependence. Figure 15 is a similar plot for three 
levels of £. The scatter inherent in the £ data is nicely illustrated 
by removing the height dependence in a .U^ time history (Fig. 15) of buoy 
data. Note that the single level scatter is much smaller than the level- 
to-level disagreements. The average profile estimate (based on Eq. 19) 
of is about .12 m/sec whereas the average turbulence estimate is about 
.10 m/sec. This is, in fact, good agreement and illustrates the superior- 
ity of the buoy system for gradient measurements. The compatability of 
gradient and turbulence measurements can be examined in calculations of 
normalized heat flux, F^, using the profile method 

Fh = - 06/3Z) (20) 

and the turbulence method 

where and F(R^) are functions of stability. These methods have been 
favorably compared (Fairall, et al., 1978) for fog events during a marine 
fog study called CEWCOM-76 (Fig. 16). 



C. Summary 

It is apparent that this line of research is permeated with 
instrumental difficulties. The atmospheric gradients over the ocean are 
small enough so that even state-of-the-art mean measurements give fairly 
poor accuracy. The turbulence measurements are inherently well above 
noise levels but the sensors are subject to an array of environmental 
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problems; salt and water films, platform motion, aerodynamic flow dis- 
tortion by the ship, radio and radar interference, and even corrosion. 
Despite these problems, the system described in this paper has yielded 
valuable and unique data about the open ocean boundary layer. Those 
theoretical discrepancies that still exist will hopefully be resolved 
through the continual evolution and improvement of instrumentation and 
measurement techniques. 
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Figure Captions 



1. Naval Postgraduate School Department of Oceanography's R/V Acania, 

2. Buoy system for near surface measurements. 

3. Block diagram of mean system. 

4. Schematic of aspirated shelter and sensors. H: humidity sensor / 

Q: quartz thermometer, 0: quartz oscillator, F: exhaust fan, 

R: radiation shield. 

5. Block diagram of fluctuation system. 

6. Wind tunnel calibration of mean wind sensors (cups), sensor error vs. 
wind speed . 

7. Lyman-alpha calibration, output voltage vs. water vapor pressure - 

8. Example profiles of temperature, humidity and wind speed. 

9. Example power spectra of fluctuation parameters vs. frequency, 

the low frequency roll-off is instrumental ac coupling, the high frequency 
roll-off represents the beginning of the dissipation range: the slanted 
lines represent Kolmogorov -5/3 exponential dependence a) temperature , 
b) humidity, and c) velocity. 

10. Comparison of values from the spectral method and the RMS method. 

11. Route of R/V Acania to and from SNI. 

12. Positions and track of R/V Acania in vicinity of SNI. 

13. Anchorage locations in Monterey Bay. 

14. Turbulence parameters vs. local time during CEWCOM-76 for Z = 4.2 m 
(open circles) , Z = 6.6 m (solid circles) , Z = 11.3 m (triangles) , 
and Z = 17.7 m (x's): a) C^^ and b) £. 

15. Friction velocity (U^) vs. local time comparing three turbulence levels 
and the gradient value . 

16. Calculation of sensible heat flux, W = C Fj^ where C is the heat capacity 
of air, using the mean gradients (open circles) and the turbulence 
values (x's) . 
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APPENDIX B 

MOS STABILITY FUNCTIONS 



The forms of the mean gradient functions (Businger et al., 1971). 



- (1 - 


5<0 


4.^(5) = (1 + 4.7c) 


OO 


(^^(C) = (1 - 19C)‘''/^ 


5<0 


<t>j(c) = (1 + 6.4c) 


OO 


The mean profile function 


2 

i>U) = l-i,(c) - 31n ^(c) + 21n (l.^-^ii-0.) ^ ? tan"'' ^(c) - ^/2 + In (^^^^^) 


The dimensionless velocity dissipation function 


(Wyngaard and Cote, 1971) 


E^(C) = (1 + .51c^/^)^/2 


OO 


E^j(c) = (1 + 2.5c 2/3)372 


OO 


The dimensionless temperature structure function 


parameter (Wyngaard et al, 


1971). 


f.^(c) = 4.9(1 - 7C 


C<0 


f.p(c) = 4.9(1 + 2.4C 2/3) 


OO 



B-1 



APPENDIX C 
DATA 



C-1 



Table I-C. EM Propagation data (bulk calculations). 



Z* duct height (m) 

Fm Minimum Frequency tripped (GHz) 

N Refractivity (N Units) 

dN/dZ Refractivity gradient (tn~^) 

C|^ Refractive index structure function (m""^' ) 



C-2 



Date 



Z* 



Frr. 



d;s/o, 



Cr"2 



-j— 1 2 /2r4 F4tT9 i- 5 - 7-6 3 3 1 . 3 ■rt-;-2-4-4 i— 3 9F — Y> 

2 2/23 1500 17.4 5.0 332.7 -0.24'^ 1.04E-12 

3 2/23 1530 19.2 4.3 332.3 -0.254 1.160-12 

r— 4 2-/2-S ^l:8-r7 ^4-75 33 3 ; 0 =-0t2-5o 1.116-3^12 —- 

■ 5 2/23 1630 17.0 5.1 335.1 -0.238 8.93S-13 

’ 6 2/28 1700 16.8 5.2 334.6 -0.235 9.06E-13 

—7 2-/2 S ^t-7-3-0 3r7-r6 -5“rl 3 3 1 ; 7 •-=-^>t 34-0 0-r43'6^13- — 

8 2/20 1300 13.8 7.0 335.4 -0.207 7.27E-13 

10 3/01 1330 17.3 5.0 324.5 -0.277 2.04E-12 

r-ti ^3 / 0-i 1 90 0 1-6 t 3- 533 334-r5 ^0 t2T1 2r32 E -12 — 

12 3/01 1930 17.4 4.9 322.4 -0.232 2.19F-12 

13 3/01 2000 17.1 5.1 322.9 -0.273 2.15E-12 

^3/t/i -20-30 ^Jr6-r7 — “ "S .3 334-;'4 1 tO"7 E - 12 

15 3/01 2100 17.5 4.9 324.3 -0.231 2.03E-12 

16 3/01 2130 17.5 4.9 323.3 -0.232 2.15E-12 

"It? 3 -/O-I— — -2200 12t-3 5t0 322-.^} ~tyr2T9 2'rl2E:=‘12 

18 3/01 2230 16.9 5.2 324.4 -0.273 2.01E-12 

119 3/01 2300 17.3 5.0 324.6 -0.277 2.01E-12 

-i-20 3/0-1 -23-30 ^10-; 0 3 t2 325-76 ^ 0 7 26 9 IrrSSE-^ 

21 3/01 2400 15.0 5.8 329.0 -0.246 1.60E-12 

22 3/02 30 15.3 6.0 329.9 -0.241 1.42E-12 

•"23 -3/’02 Itro rSTTl 6 t 1 330 .9 =^J”.236 ' rr32E- 13 — 

24 3/02 130 14.2 5.3 332.9 -G.220 l.llE-12 

125 3/02 200 13.3 7.0 333.9 -0.214 1. 023-12 

•^23 -3-/02 230 r 3 T 2 rro 333-; o- ^3 .-393 3 .~ 33 ~-r 3 ■ ■ 

27 3/02 300 12.8 7.3 3 35.9 -0. 193 3.22E-13 

28 3/02 330 12.6 3.1 336.1 -0.195 7.98E-13 

'•h2-9 -2/-0-2— 4 0-0— 12 t3~ 3t4 323^75 "^32290 7'. 5562-13 

130 3/02 430 12.5 3.2 336.5 -0.193 7.59E-13 

21 3/02 500 11.7 9.0 337.3 -0.132 5.55E-13 

'22 3/0-2 5-30 tl .“1 9t7 -333t3 =197-373- 3-73252=33 

23 3/02 500 10.5 10.5 340.1 -0.155 4.97S-13 

'34 3/02 630 3.9 13.5 342.2 -0.142 3.53E-13 

-32— 3/02 -7x90 3 t2 23t 6 34-4-74 ^-33:3 27-r5Er-r3 

36 3/02 730 5.0 32.7 346.4 -0.037 1. 368-13 

37 3/02 300 4.5 35.4 347.0 -0.0G4 l.lSE-13 

23 3-/32 3-33 3t2 52-2 -342-73 =0-72-72 Trl'- E= 1 4 

39 3/02 900 3.3 59.5 349.0 -0.070 5.45E-14 

40 3/02 930 3.7 50.1 349.1 -0.G75 6.43E-U 

-43: 3/22 122^3 472 42r2 34272 ^tTrOT/t! o-r3-3o-12 - 

42 3/02 1130 4.4 38.4 349.3 -0.087 7.29S-14 

43 3/02 1200 5.0 32.3 348.1 -0.093 9.60E-14 

-22 3-/3 1 1-220 772 42t 2 342-.-2 -07 0 7 2 5-r27rr--14 — 

45 3/02 1300 2.3 75.7 350.4 -0.065 3.675-1* 

46 3/02 1330 2.8 75.3 350.2 -0.064 4.23E-14 

-4-7 3/0 2- 3420 311-73 351.0 =27237 2-7922=12 

48 3/02 1430 2.6 86.7 350.2 -0.062 3.89E-14 

49 3/02 1500 3.5 54.7 349.1 -0.073 6.32E-14 



# Date GMT Z* Fai N dl'J/dZ Cn''2 



—5 0 3/0-2 153 0 5r2 54-. -0- 34 0-^8 ^065 4-r95B- 14 

51 3/02 1500 2.7 81.3 350.3 -0.064 3.56E-14 

52 3/02 1530 3.7 45.7 343.1 -0.073 8.10E-14 

^3/02 — 170-0 -2-i-5 01 t5 34-5-.-7 OtOS-? ^3-55B-14 — 

54 3/02 1330 3.0 58.1 350.5 -0.071 2.99E-14 

55 3/02 1900 3.5 51.9 349.4 -0.075 5.15S-14 

— 55 3/-0-2 193 0- 2-.-0 i-tl5-r5 35-i-rO 1 -. 8-75-1. 4 

57 3/02 2000 3.1 66.8 350.1 -0.070 3.27E-14 

53 3/02 2030 2.7 79.3 350.7 -0.G58 2.21E-14 

—59 3/0-2 -2-iOt) 2-r-T 51-r9 550x9 O-rOrT 2-.05 1-1 4 

60 3/02 2130 2.7 81.9 350.9 -0.057 2.23E-11 

61 3/02 2200 2.3 78.6 350.7 -0.053 2.33E-14 

-5-2 3-/92 22-30 3x9 6-2x5 34-0.-9 -8xO-7-3 ^3-.-4-9-E— H 

63 3/02 2300 3.1 54.4 350.1 -0.072 3.34E-14 

64 3/02 2330 4.5 33.0 348.5 -0.035 5.93E-14 

-05- 3tO- 2 240 G- 6x0 24-r-2 34-6x9 Ox-ldO -t-r-2-4-F--t-3 

55 3/03 30 6.3 22.7 345.4 -0.111 1.37E-13 

57 3/03 100 5.7 20.9 346.2 -0.115 1.54E-13 

-0-8 -3/-03 1-30 6-xS— 20x9 34-5 xO 0x11-7 0x4 4-5^13 

69 3/03 200 3.1 15.5 344.5 -0.133 2.25E-13 

70 3/03 230 12.6 7.3 333.6 -0.194 6.22E-13 

— 7-1 3 740 3 300— 15x9 6-x9 3-30r4 —0x006 070413-02 

72 3/03 330 13.9 6.9 326.9 -0.220 1.65E-12 

73 3/03 400 13.5 7.3 328.5 -0.214 1.41E-12 

—74 3-/ 0-3 4-3 0 2 0.-0—- ^3-r9 3-2-3-r9 -0x30-0 — — -2-r3-2E -iO 

75 3/03 500 18.0 4.7 328.2 -0.275 1.71E-12 

76 3/03 530 13.3 4.6 325.9 -0.235 1.99E-12 

—77 0-/03 600 13-r2 4x6 305x-5 -0x005 ^2 .-0-3E-i2 

78 3/03 530 20.3 3.3 319.6 -0.333 2.94E-12 

79 3/03 700 20.4 3.5 318,6 -0,329 3.02E-12 

—80 -0-/03 730 19-r3 4x0 300x9 OhrOOS -2x6-3-5 - 12 

81 3/03 800 17.9 4.7 323.7 -0.235 2.21E-12 

82 3/03 830 20.6 3.8 317.9 -0.334 3.13E-12 

— S-3 3/0-3 960 ^20;-9 8-78 300x0 OxlOO 3x-10B-i2 

84 3/03 930 19,1 4,3 319.7 -0.309 2.75E-12 

85 3/03 1000 19.3 4.2 321.4 -0.311 2.59E-12 

—36 3 /0-3- ^1000 t6x-l 5x6 300x1 '0x050 2x10-0-12 

87 3/03 1100 13.5 4.5 320.3 -0.301 2.53E-12 

S3 3/03 1300 15.2 5,5 325.2 -0.258 1.S6E-12 

—09 — 3/-0-3 -103 0 ^F7-x4 5~r 0 - 3 2 -4- . 8 -Ox-206- -2x1-00-1-2 

90 3/03 1500 20.2 4.0 322.4 -0.317 2.52E-12 

91 3/03 1630 19,5 4.2 323.9 -0.305 2.31S-12 

—90 -3A73 1700 OOtO 3x0 002x0 - 0. 318 2 . 56 -E- 1 2 

93 3/03 1830 14.2 5 .8 3 30 .3 -0.221 1.37E-12 

94 3/03 1900 15.7 5.8 323.9 -0.246 1.56E-12 

—96 3 1-3-30 ^16x0 5x0 OOOxO -0x069 1-.-566-10 

96 3/03 2000 14.4 5.6 328.5 -0.225 1.51E-12 

97 3/03 2030 14.0 5.9 329.7 -0.220 1.34^-12 



6-4 



Date G^'.r 



2* Fr l-i dN/dZ Cr'‘2 



-•9^ -3t^> 3 2 IDO -}-6-r ‘T 5-7-6 3 - 2 " 5 .0 ~t)— 2D2 1 12 

j93 3/03 2130 13.7 7.1 329.8 -0.216 1.32E-12 

00 3/03 2200 12.3 8.3 332.2 -G.193 1.05F-12 

Di- -3-/Q3~ 2 - 23 -e t6-r7 5 -t3 3-2-3T9 ^t!v2tjtr 2r0l2-12 

02 3/03 2300 17.0 5.1 323.3 -0.272 2.04E-12 

03 3/03 2330 15.2 5.1 327.7 -0.240 1.56--12 

-04 3/0-3 2400 17rO 4-7-7 3-22-r3 --Ov-20-7- 2.-2aD-12 

05 3/04 30 18.3 4.6 320.2 -C.237 2.57F-12 

06 3/04 100 16.4 5.4 323.5 -0.264 2.03E-12 

0-7 3/04 30 irhr5 4et 9 O-O-hTl ^^-0-r28O 2 -36 E^!r2 

,03 3/04 200 17.2 5.0 320.3 -0.230 2.37F-12 

00 3/04 230 13.3 4.5 319.5 -0.296 2.60S-12 

■^445 3/04 ^400 3r7-rO ttI 30:0^7 --Orr2-70 2*70 0^>-12 

.11 3/04 430 15.6 5.3 319.5 -0.270 2.35S-12 

.12 3/04 500 16.9 5.2 318.9 -0.275 2.44E-12 

*tl-3 3/0 4 5-3 0 ^t6T5 5 tt 3 3 tOrO ^Ot 2-72 2-.-3 9 I2r 

.14 3/04 500 16.4 5.4 318.9 -0.268 2.34E-12 

.15 3/04 530 17.3 5.0 318.5 -0.282 2.55E-12 

-jOrO — —3/04— 700 l-7vT 4-rO 3-1:074 =0 t201 2tOOE-^ 12 — 

117 3/04 730 17.3 5.0 316.1 -0.234 2.75E-12 

118 3/04 800 '16.7 5.3 318.3 -0.273 2.47E-12 

-rl-9 37434 3-30 5 t 3 3 TOt 4 ~(rrZ 7 X) 270-4F=^r2 — 

120 3/04 900 15.9 5.2 319.3 -0.276 2.44E-12 

L21 3/04 930 17.1 5.1 319.7 -0.279 2.43E-12 

-t-2-2 2/0-4 1 0 00 33779 6-79 32 - ^ . 9 =337223 Iv 06 E=^12 

123 3/04 1030 15.9 5.7 321.8 -0.256 2.13E-12 

124 3/04 1100 17.2 5.0 319.7 -0.281 2.50E-12 

»325 ^3-/0-4 3330 377“4 4T3 321-70 =0t28-2— —2^350=32 

126 3/04 1330 17.0 5.1 319.2 -0.273 2.50E-12 

127 3/04 1400 18.4 4.6 321.4 -0.297 2.50E-12 

'i-23 3-/34 14-30 377-5 4-79 33072 =0 t 23O 2707r-=l 2 

129 3/04 1600 13.5 4.5 317.6 -0.303 2.35E-12 

130 3/04 1800 16.7 5.3 313.9 -0.272 2.42E-12 

1333 37^0-4 3 '303 38772 4776 -335-74 =0-7200 ' 2T '"Or - 32 

132 3/04 1930 13.2 4.7 315.1 -0.299 2.95E-12 

133 3/04 2000 13.1 4.7 315.1 -0.299 2.54F-12 

^33-^ ^3/04 2-330 ^33t 3 trS 337 4 t 2 =077-30-3 3rtr4 — -]2 

135 3/04 2100 17.9 4.7 315.1 -0.296 2.91E-12 

135 3/04 2130 19.0 4.3 312.4 -0.316 3.34F-12 

f±33 3/0-4 -2200 ^1:272 472 3O0 t 6 =---07222 ^37-6-tF=-l-2 

133 3/04 2230 13.4 4.5 312.2 -0.306 3.2^E-12 

139 3/04 2300 13.4 4.6 312.3 -0.305 3.22E-12 

3t 4-0 ^3/24 23-30 l"7-r6 4-rO =0 t2O0 2t24- — 22 

141 3/05 100 16.3 5.5 315.5 -0.253 2.52E-12 

142 3/05 130 15.8 5.3 316.7 -0.258 2.3^7-12 

± 4-3 ^3/05 2'00 2575 5-7 O 32275 --27 2 5 -4 27729 et=-12 

144 3/05 300 15.2 5.0 315.1 -0.249 2.26E-12 

145 3/05 330 15.1 6.1 315.2 -0.247 2.24E-12 



Date 



Fn 



cl K/a ; 



Cn' 



3-/ 0^ 4 00 h5-r-l -O-r-2 30:4^4 ^—2^6 2 5 F -i ^ — 

' 147 3/05 430 15.2 5.0 313.0 -0.250 2.34F-12 

143 3/05 500 15.4 6.0 312.5 -0.253 2.405-12 

^ 9 2/Q 5 5 3 0 i-4-r0 -5^5 3 07-2-27 1 . 9 3 '1 - 1 2 

150 3/05 500 14.1 6.8 311.8 -0.230 2.00E-12 

151 3/05 630 13.4 7.3 309.8 -0.219 1.79E-12 

1-5 2 -3/-05 700 16-rO 7-;7 309-v9 0.-2 tl -1-725-12 - 

153 3/05 730 12. S 7.9 308.3 -0.207 1.64E-12 

154 3/05 800 13.0 7.7 306.8 -0.211 1.69E-12 

—1-35 ^3-/05 8-30 -i-M 8-r3- 30-0^5 -0-.-i90 t.-^7E--i2 

156 3/05 900 12.1 8.6 303.6 -0.195 1.33E-12 

157 3/05 930 12.1 8.6 307.1 -0.195 1.40S-12 

JrO-l 3-/0-5 t03O 3 t2^0 3-t 6 300-r^^ -0-.-3r05 ^]r.-305-i-2 

162 3/05 2000 23.5 3.2 311.5 -0.384 4.59E-12 

153 3/05 2030 20.5 3.9 312.2 -0.340 3.90E-12 

: ^104 3-/0-5 2-V0O 20-7-3 3-70 31:0-r7 4-.-0->E'>0:2 

165 3/05 2130 21.3 3.7 310.1 -0.355 4.33S-12 

155 3/05 2200 21.0 3.7 310.4 -0.351 4.26E-12 

10-7 3/0-5 -2-200 2-/t 0 3-r4 -00-07-2 -0r-300 1-2 

163 3/05 2300 20.0 4.0 312.0 -0.333 3.S9E-12 

159 3/05 2330 20.0 4 .0 312.0 -0.33 3 3.89E-12 

170 3-AO0 2400 20tO 4t 0 SAtri - 0 - .-331 

171 3/05 30 21.0 3.3 312.3 -0.341 3.67E-12 

172 3/05 100 21.3 3.7 311.1 -0.349 3.36E-12 

^173 3-A30 1-08 ^107-6 ■ 4 7'1 04070 07-004 07-06-0- 10 

174 3/06 200 22.4 3.4 310.1 -0.355 4.20E-12 

175 3/OG 2 30 21.5 3 . 6 308.4 -0.357 <i.30E-12 

^t-70 0A>6 300 03r . 5 3 7 5 3 8-0 t 0 -0 . 356 4-7-26-E -ti 

177 3/06 330 21.1 3.7 303.5 -C.351 4.38E-12 

173 3/05 400 21.0 3.7 309.7 -0.351 4. 540-12 

77-9 -3/0 6 400 0070 3-r7 300 7 ^ -9 - . - 0 - 5 - 1 4-r600-00— — 

180 3/05 530 20.2 4.0 311.7 -0.336 4.14E-12 

131 3/06 500 20.4 3.9 313.2 -0.336 3.87E-12 

03-2 3/-05 630 3r07r2 070 00070 - 07 - 317 070-4f— t2 

133 3/05 730 17.3 • 4.3 317.2 -C.239 2.S0E-12 

196 3/07 730 12.6 S.l 320.7 -0.201 1.415-12 

^777 -3A37 000 72-r0 3 -t0 0 01 - -4 - 0—193 07000—10 — 

200 3/07 930 15.5 5.3 320.7 -0.253 1.93E-12 

201 3/07 1000 14. S 5.3 320.3 -0.237 1.59E-12 

00-2 3-/07 10-30 ^107-6 4-r9 02070 -0^076 0r00E-00 

203 3/07 1100 16.8 5.2 324.6 -0.250 1.65E-12 

204 3/07 1130 17.0 5.1 324.0 -0.268 1.91E-12 

20-5 -0/0-7- — 1200 76-. 4 5-r4 3 - 27 - 7 - 7 -0-202 17-05&-02 

206 3/07 1300 14.9 6.2 329.5 -0.232 1.39E-12 

207 3/07 1500 8.3 14.9 327.5 -0.131 4.98E-13 

, 2-0 8 3/0-7 7000 0^0 10r0 306-r8 -07-000 4. 58 E - H 

.! 209 3/07 1800 7.4 17.7 332.5 -0.117 3.68E-13 

211 3/07 1900 8.0 16.0 332.5 -0.125 4.01E-13 



Date GMT Z* Fm N dN/dZ Cn"2 



•L-2 3-/07 1930 8-r6 i^^-2 3-29r9 O^i-35 5rl2E-l-3 - 

,L3 3/07 2000 8.6 14.3 329.6 -0.135 5.14E-13 

:14 3/07 2030 9. 5 12.4 326.4 -0.149 6.41E-13 

iir5 -3-/0-7 2100 lir-2 9%-7 ^3-2-3r-2 9.59E-13 

:i6 3/07 2130 9.3 12.7 324.0 -0.147 6 .09E-13 

:l7 3/07 2200 10.2 11.0 323.5 -0.161 7.55E-13 

it8 3/07 2-239 1-0—2 idrrO — — -32-2“,-8 “'Thr'1-6-2 ^tT-OE-IO — 

19 3/08 130 14.3 6.7 322.8 -0.230 1.64E-12 

9.0 3/08 730 12.9 7.7 317.9 -0.209 1.66E-12 

|23r 3/00 300 l-6-r8 5~2 3'lrOT4 -0'7 2715 2t7 0E-1-2 

9.2 3/03 830 16.9 5.2 318.3 -0.276 2.67E-12 

p 3/03 900 17.7 4.8 3 17.3 -0.291 2.97r,-12 

B-4 3-/0 3 9 3u -i-3-nl rr± 3 ir 7 - ;6 -07^3 1 2 

25 3/08 1100 14.5 6.5 313.2 -0.236 2.15E-12 

26 3/03 1130 1G.9 5.2 315.5 -0.253 1.52E-12 

^7 3/0 3 1-2-3;2 2±r^ /tS 3ir3r-3 -^3^74 dr-4-± E^12' 

j23 3/03 1400 20.3 3 .9 315.1 -0.336 3.95C-12 

29 3/08 1430 13. S 4.4 315.4 -0.314 3.62E-12 

!30 ^3/9-s iii-OO illri ‘tTT 3ir6T"2 =1^r30tr SttIOT-I-I 

i31 3/03 1530 14.3 6.7 314.4 -0 . 234 2.17F-12 

32 3/OS 1700 13.7 7.1 315.3 -G.224 1.940-12 

^3 -3:9- 30 13 .- 5 Tr^ 316.4 ^^ ^ O t 3 22 -1 ^iTtrTlT^Irr 

34 3/08 2000 18.1 4.7 317.0 -0.300 3.1SE-12 

35 3/03 2030 11.7 9.0 318.1 -0.133 1.29E-12 

3^ y/ITc 2r2tyu tlr-? TrC 31r7vr -0-rl93-8 ^rr2^E--="12 

37 3/03 2230 11.7 3.9 313.6 -0.139 1.31E-12 

30 3/03 2300 11.7 9.0 316.- -0,137 1.29E-12 

41 379x9 tO'09 10 . 5 -“It) 75 31^3 "-BUTTS'S IVCnOET-' 12 

42 3/03 1030 14.6 5.5 320.4 -0.237 2.05E-12 

43 3/09 1200 9.7 12.0 323.0 -0.152 7.81E-13 

44 rtro 975 -12-73 mr9 ^tTTTSO T.23r=d.-3 

45 3/09 1300 10.0 11.4 321.5 -G.158 3.33E-13 

46 3/09 1400 18.0 4.7 322.3 -0.289 2.47E-12 

•47 3-/0^ "1420 rrro ^474 321175 — -tr73-cr3 rrs orr~ 12 

48 3/09 1500 19.5 4.2 321.7 -0.310 2.790-12 

49 3/09 1530 13.6 7.1 321.9 -0.220 1.710-12 

513 37999 lo20 ^r8T4 475 322 ~ .2 =47245 -2-.“r7 0 - 12“ 

51 3/09 1630 13.3 4.4 323.7 -0.299 2.53E-12 

52 3/09 1700 18.5 4.5 324.4 -0.203 2.470-12 

33 3/09 lirtK) rW 5-75 325 t3 ^47253 2". 104-12 

54 3/09 1330 12.9 7.3 326.0 -0.207 1.500-12 

55 3/09 1900 11.6 9.1 328.2 -0.134 1.12E-12 

46 3/ 09“ 20ryd 12rr3 749 329-74 =94904 14-32 9-=12 

57 3/09 2100 12.7 3.0 327.9 -0.202 1.31S-12 

58 3/09 2130 12.3 7.3 329.2 -0.204 1.23r-i2 

59 3/0-9 2200 1977 749 326.9 =44229: 1. 51 0 = 19 - 

60 3/03 2230 14.6 5.4 325.7 -0.235 1.730-12 

61 3/09 2300 13.4 7.3 327.9 -0.215 1.43E-12 
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-26^ d-/l-0 S-36- i ih-.a ^-7^ -y-.~y3^ •]tt-2'3-5~ 12- 

263 3/10 700 15.1 • 5 .2 331 .2 -0.23'’- 1.46E-12 

254 3/10 730 15.3 5.0 331.2 -C.2’2 1.49E-12 

2 ^ 5 — ^ — i'/itj ^30 IrdT^- 5~7 33d". I— ^ “x>. d-'vd -jT.'dd'S'- 12 

265 3/10 900 15.3 5.7 320.5 -0.252 1.55S-12 

257 3/10 1100 15.5 5.5 335.0 -0.240 1.23E-12 

^2^3 3-/±-e 1 i-5-r-l 5-rl 3-3-5— 1 5 l-.-2-i-&-“-t2 

25b* 3/10 1200 15.4 5 .0 335.7 -0.237 1.195-12 

270 3/10 1230 14.6 5.5 337.5 -0.223 l.OOE-12 
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272 3/10 1500 12.1 8.6 311. 3 -0.187 6.73E-13 

273 3/10 1530 11.5 9.3 343.1 -0.177 5.45E-13 

2^4 3-/-10 ye 00 3t 0 4-e-r9 35-t-rO -O-rO-7-2 8-7-5-7-&-1 t4 

275 3/10 1830 4.5 37.3 350.9 -0.082 9.92E-14 

276 3/10 1900 4.3 40.2 351.5 -0.030 8.32E-14 

2-7-7 3/tO 103 0 4 t-2 y±~r2 3-3±-r7 ^OarOilO- 77-74 14 

278 3/10 2000 4.3 41.0 351.8 -0.033 6.85E-14 

279 3/10 2030 3.7 50.9 352.5 -0.075 5.25E-14 

2BQ 3-/ 19- ^21^9 ^3-r3 69-ri 35 2. - 7 ^0-r9--M 3-r99-E---14 — 

281 3/10 2130 3.7 51.2 352.5 -0.075 5.06E-14 

282 3/11 700 7.1 18.8 343.0 -0.114 3.28E-13 

20-3 3-/-1 1 7343 TtO 18-r3 3-4^-r3 -0 -7- 11 8 3-r40e---13 

284 3/11 600 7.8 16.4 341.3 -0.124 4.09E-13 

285 3/11 330 3.0 15.3 341.4 -0.127 4.25E-13 

-200- 2A-1 70 0 8-7-3 15-r2- — 3-4-17-6 -O-T-l-lt ^30-E3-ir3 

237 3/11 930 11.0 9.3 336.4 -0.174 3.31E-13 

283 3/11 1000 13.7 7.1 331.0 -0.219 1.37E-12 

23-9 3/-11 1 e-33 i^rrd 3-rr 334-rr9 --d-rirOa 1^- 12 

290 3/11 1100 11.3 9.5 335.3 • -0.179 9.05E-13 
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233 3-/H 1-233 ^ti-r3 TrO 333- r l - 0 t 3 - 31 0 - r2 2 E-13 
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294 3/11 1330 16.5 5.4 324.2 -0.259 2.21E-12 

2 b* 5 ■ — 3-/-H:— — 14-O-0 ^ 15 - 7-5 5 - 7-9 323-t3 “0-.-2 53 2-7333— '±2 
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--^3-7 3/33 2333 177-5 4-t 5 32373 -- 3 r-2: 3 3 2 t€-23=3-2 

308 3/11 2100 15.2 6 .1 326.7 -0.245 1.8'^E-12 

309 3/11 2130 15.5 5.9 324.0 -0.252 2.10E-12 
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Table II-C. MOS scaling parameters and EM propagation results. Compari- 
son of the bulk method (first line) with turbulence U* and 
profile T*, q* (second line). 



Z/L Monin-Obukhov stability parameter, 5 

U* Velocity scaling parameter (m/sec) 

T* Potential temperature scaling parameter (°C) 
q* Water vapor mixing ratio scaling parameter (gm/kg) 
Z* Duct height (m) 

Fm Minimum frequency trapped (GHz) 
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Table IV-C. Meteorological data. The data is given by level number 
(ie., U2 is the wind speed at level 2), 

U Relative wind speed (m/sec) 

PHI Relative wind direction (degrees) 

TS Sea surface temperature (°C) 

T Air temperature (°C) 

H Relative humidity (%) 
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